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Abstract 


The tubular SOFC generator CHP-100, built by Siemens Power Generation (SPG) Stationary Fuel Cells (SFC), is running at the Gas Turbine 
Technologies (GTT) in Torino (Italy), in the framework of the EOS Project. The nominal load of the generator ensures a produced electric power 
of around 105 kW. ac and around 60 kW, of thermal power at 250°C to be used for the custom tailored HVAC system. 

Several experimental sessions have been scheduled on the generator; the aim is to characterize the operation through the analysis of some global 
performance index and the detailed control of the operation of the different bundles of the whole stack. 

All the scheduled tests have been performed by applying the methodology of design of experiment; the main obtained results show the effect of 
the change of the analysed operating factors in terms of distribution of voltage and temperature over the stack. 

Fuel consumption tests give information about the sensitivity of the voltage and temperature distribution along the single bundles. 


On the other hand, since the generator is an air cooled system, the results of the tests on the air stoichs have been used to analyze the generator 
thermal management (temperature distribution and profiles) and its effect on the polarization. 

The sensitivity analysis of the local voltage to the overall fuel consumption modifications can be used as a powerful procedure to deduce the 
local distribution of fuel utilization (FU) along the single bundles: in fact, through a model obtained by deriving the polarization curve respect to 
FU, it is possible to link the distribution of voltage sensitivities to FC to the distribution of the local FU. 

The FU distribution will be shown as non-uniform, and this affects the local voltage and temperatures, causing a high warming effect in some 
rows of the generator. Therefore, a discussion around the effectiveness of the thermal regulation made by the air stoichs, in order to reduce the 
non-uniform distribution of temperature and the overheating (increasing therefore the voltage behavior along the generator) has been performed. It 
is demonstrated that the utilization of one air plenum is not effective in the thermal regulation of the whole generator, in particular in the reduction 
of the temperature gradients linked to the non-uniform fuel distribution. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Gas Turbine Technologies (GTT) and Politecnico di Torino, 
both located in Torino (Italy), have installed a SOFC laboratory 
in order to analyze the operation, in cogenerative configuration, 
of the CHP-100 SOFC Field Unit built by Siemens Power Gen- 
eration (SPG) Stationary Fuel Cells (SFC). The generator was 
commissioned on June 19, 2005, and to date has shown the 
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record availability of 99.8%, showing more than 30,000 operat- 
ing hours at this scale. The demonstration of this generator is the 
first phase of the EOS project (2004-2009), that has the follow- 
ing aims: (i) to establish a research group that, joining university 
and industry, will allow the growth of the scientific knowledge 
about SOFC technologies; (ii) to provide GTT with installation 
and operational know-how with SOFC; (iii) to extend life testing 
of SOFCs of varying pedigrees; (iv) to perform experimental 
and modeling activity to characterize the generator behavior. 
At present, approximately 100 kW, of ac electrical energy are 
delivered to the grid, and 55kW, of thermal energy provides 
heating and cooling for the facility’s buildings upon integration 
with a district HVAC. Politecnico di Torino has developed sev- 
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Nomenclature 

BoP balance of plant 

CHP cogeneration heat and power 

Di-j fundamental binary diffusivity in the electrode 
(cm? s7!) 

Daeft) effective binary diffusivity in the anode layer 
(cm? s7!) 

Deef effective binary diffusivity in the cathode layer 
(cm? s7!) 

Eo open circuit voltage (V) 

F Faraday number (C mol`!) 

FC fuel consumption (%) 

FU fuel utilization (%) 

HC high current 

las anode limiting current (A cm~?) 

leg cathode limiting current (A cm7?) 

ic cell current density (A cm7?) 

I generator current (A) 

L cell length (m) 

LC low current 

p cell pressure (Pa) 

Pity hydrogen pressure at the anode/electrolyte 
interface (Pa) 

PR, hydrogen pressure at the anode bulk (Pa) 

PO, oxygen pressure at the cathode/electrolyte 
interface (Pa) 

Pb, oxygen pressure at the cathode bulk (Pa) 

R universal gas constant (J mol! K7!) 

SOFC solid oxide fuel cell 

ta thickness of the anode layer (cm) 

te thickness of the cathode layer (cm) 

T temperature (K) 

Togn setpoint temperature (K) 

Tair air pre-heating temperature (K) 

Ve cell terminal voltage (V) 

Vaiff diffusion overpotential (V) 

Vnernst Nernst potential (V) 

VNernst position averaged Nernst potential (V) 


VNernst,HC average Nernst potential, high current 


approximation (V) 


Vnernst,LC average Nernst potential, low current 


Xi 
Vi 
yip 
Yep o 
Y; 


approximation (V) 

axial position along the tubular cell (m) 
regressor variables 

molar fraction of reactants at the electrodes’ bulk 
hydrogen molar fraction at fuel inlet condition 


water molar fraction at fuel inlet condition 
dependent variable 


Greek letters 


Bi 
E 
Nact,alc 


regression coefficient 
electrode porosity 
activation overpotential at anode and cathode (V) 


Nconc,ae concentration overpotential at anode and cathode 


(V) 
Nohm ohmic overpotential (V) 
Aair process air excess 
T electrode tortuosity 


eral activities, like the design of the SOFC primary generator 
thermal plant [1], the modeling of the generator and of the BoP 
[2-5], the safety design and analysis [6]. 

In particular, Politecnico has developed several experimental 
test and analysis of the experimental results. The experimental 
sessions have been designed in order to investigate the effect 
of four control factors at the nominal load operation: the fuel 
consumption FC, the air stoichs Aair, the setpoint temperature 
Toen, and the air pre-heating temperature Tair. The methodology 
of the Design of Experiments has been used in the planning of 
the experimental session (factorial design) [7,8]. In [9] a model 
of the operation of the SOFC CHP 100 has been developed 
using a OD approach and validated through a first session of 
experimental tests. In [10,11] the effect of the setpoint genera- 
tor temperature and fuel utilization factors on several dependent 
variables (i.e. dc and ac electric power, recovered heat, electric 
and global efficiency, efficiency of the pre-reforming process) is 
analyzed in form of screening tests and the process responses are 
treated in form of response surface plots. In [12] the regression 
models have been used in constrained optimization procedures, 
to have an overview of the optimal operation points which max- 
imize different objective functions (ac electric power, recovered 
heat, etc.). 

The experimental tests showed that the most significant fac- 
tors affecting the generator behavior are the fuel consumption 
FC and the air stoichs A,i,. In this paper results are presented 
concerning a deeper analysis of the sensitivity of voltage and 
temperature distribution to modification of fuel consumption 
and air stoichs. The aim is the deduction of the local distribution 
of fuel utilization (FU) along the single bundles, its drawbacks 
on temperature distribution (possible overheating), and a dis- 
cussion concerning the temperature regulation through the air 
stoichs which outline some problems. 

In [13] the sensitivity of the measured terminal voltage of 
a tubular solid oxide fuel cell to the fuel utilization is used to 
provide information about ‘leaks’ through the cell. In [14] an 
experimental analysis on a short stack is presented and the pres- 
ence of localized overheating is addressed to low fuel availability 
on the anode side. 


2. Experimental 
2.1. Description of the plant 


The CHP-100kW,. SOFC Field Unit (SPG SFC) utilizes 
the commercial prototype air electrode supported cells (22 mm 
diameter, 150 cm active length, 834 cm? active area) and in-stack 
reformers. The generator is fed with natural gas from the grid. 
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Fig. 1. Picture of the SOFC CHP test site in GTT (Torino, Italy). 


In Fig. 1 the picture representing the CHP 100 test site in GTT 
(Torino, Italy) is shown. 

The upper level of system hierarchy after the single cell is 
the cell bundle, which consists of a 24-cell array arranged as 8 
cells in electrical series by 3 cells in electrical parallel. Four cell 
bundles are connected in series to form a bundle row, and 12 
bundle rows are aligned side by side, interconnected in serpen- 
tine fashion with an in-stack reformer between each bundle row 
(for a total of 1152 single cells). In the following of the paper we 
refer to a cell sector, which is defined like the electrical serie of 
two cell bundles. In the new rebuilt a multi-pedigree cell stack 
has been assembled; SOFC cells mainly differ for total running 
hours; one of the interest of the experimentation is to study how 
cell with different operating history work together in a stack. The 
rest of the system is composed by the balance of plant (BoP), 
with five major skids: Generator Module, Electrical Control Sys- 
tem, Fuel Supply System (FSS), Thermal Management System 
(TMS), and Heat Export System (HES) [1-6,9—12, 15-17]. 

From the electrical point of view, the module is made up of 
four parts: (a) the electrochemical generator; (b) the Power Con- 
ditioning System (PCS): an inverter, which operates the dc/ac 
conversion; (c) the SOFC auxiliaries (blowers); (d) the Fuel Cell 
Board (QFC), designed by Politecnico of Torino, which allows 
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Fig. 2. Schematic of the stack arrangement. 


one the electrical connection of the CHP 100 to the grid. From 
the thermal point of view, the exhausts from the stack, passing 
in a cross finned tubes gas-water exchanger, provide, in nom- 
inal conditions, approximately 60 kW; of thermal energy, used 
for the winter and summer conditioning (through a absorption 
refrigerator cycle water-lithium bromide fed with hot water) of 
some offices of the GTT factory. 

The schematic of the cell stack arrangement is shown in Fig. 2 
[17], outlining the position of the fuel ejectors and the power 
leads. The simplified flow schematic of the SOFC CHP-100 BoP 
is shown in Fig. 3, with the schematic of the primary generator 
structure. In Fig. 2, itis shown the position of the fuel ejectors and 
power leads. In the paper, the results will be discussed dividing 
ideally the primary generator in four main zones (North; South; 
Power Leads side and Ejectors side), and also in the central zone 
(middle between power leads and ejectors side). 

The commissioning date in GTT was June 19, 2005. The oper- 
ational data at September 2006 report a number of run hours in 
GTT of around 10,000, with a total run hours of around 30,000 
(including run hours in the Netherlands and Germany); a average 
stack temperature of 954 °C; adc generated power of 123.6 kWe, 
at 246.1 VDC and 502.3 A; a ac generated power of 113 kWe; a 
power to GTT workshop grid of 103 kW, (20% of the workshop 
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Fig. 3. Simplified flow schematic of the SOFC CHP-100 BoP, and schematic of the primary generator structure. 
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requirement); a heat generation of 60 kW; (hot water at 85 °C). 
Other characteristics are: a reliable operation also with signif- 
icant LHV changes of the natural gas; no measurable voltage 
or power degradation; a very high availability of 99.5% (actual, 
on annual basis); 4 stops (one operation error, three inverter 
failures); 24 successful operation transition to power dissipator 
due to utility (AEMD, Azienda Energetica Metropolitana Torino 
Distribuzione) or GTT grid failures; no thermal cycles; a very 
dependable automatic operation (no operator in control room); 
a remote control capability via modem; a easy maintenance 
(replacement of air filters and of desulphurization reactant). In 
1 year of operation the SOFC CHP-100 fuelled by natural gas 
avoids the production of 272 tonnes of CO2, the environmen- 
tal pollution of 1035kg of NO,, the import of 121 tonnes of 
oil equivalent with respect to a gas turbine plant of higher size 
(in the order of 100 MWe; therefore, if compared to a similar 
size gas turbine, the advantage would be even higher). Green 
Certificates are provided by the Italian Authority (GRTN) for 
electrical energy produced by a fuel cell generator (also using a 
fossil fuel): the energy production estimate in 2006 is 900 MWh 
(18 certificates). 


2.2. Description of the experimental sessions 


The scheduled tests have been carried out in the November 
2005 and March 2006, following the methodology of the Design 
of Experiments. The main aim has been the characterization 
of the operation of the single sectors of the system, pointing 
out the analysis on the behavior of the local temperatures and 
voltages. 

Four factors were investigated at the nominal load operation: 
the fuel consumption FC, the air stoichs air, the setpoint temper- 
ature Toen, and the air pre-heating temperature Tair. The most 
significant factors affecting the generator behavior are the fuel 
consumption FC, the air stoichs Agir. 

The fuel utilization (FU) factor is defined as the ratio of the 
fuel mass flow that operates the electrochemical oxidation on 
the anode surface, to the total fuel mass flow that enters the 
generator. The control variable fuel consumption (FC) factor is 
defined in the same way, but, at the numerator it also takes into 
account the fuel mass flow that is chemically oxidized on the 
anode side due to the air leakage from the cathode to the anode 
through the electrolyte layer or other sources (of the order of 
20-30 cm? min~! at 300 mA cm~? according to [13]), plus the 
fuel by-pass through the insulation package. 

The air stoichs factor (A,;,) is defined as the ratio of the total 
air mass flow which enters in the generator and the air mass flow 
operating in the electrochemical reaction on the cathode surface. 

The generator setpoint temperature (Taggy) is the higher mea- 
sured temperature by the five thermocouples placed in the central 
zone of the generator. 

The air pre-heating temperature (7};,) is the temperature of the 
process air after its heating in the thermal management system 
of the balance of plant and before entering the generator canister. 

In order to ensure a safe operation of the generator, the 
investigated experimental domain has been in the range of: 
81.75% < FC < 84.25% and 4.6 < air < 4.8, moreover, the aver- 


age current density has been kept constant at the nominal value 
of 0.2 A cm~? (500A generator current). 

The investigated dependent variables focused on the distri- 
bution in the generator of voltages and temperatures (average 
temperatures of the generator in different planes: total stack 
thermocouples, top plane thermocouples, middle plane thermo- 
couples and bottom plane thermocouples). From an electrical 
point of view it is possible to acquire the generator terminal 
voltage and the voltage of two bundles in series (sector). 

Once the setpoint was established, the data collection started 
after 60 min (to stabilize the signal data) for a period of 30 min, 
with a rate of 1 min. The single value of a variable is the mean 
value of the 30 values collected in every test length. To develop 
the study, the factorial analysis and the response surface method 
(RSM) have been applied, and first and second-order regression 
models linking the dependent variable to the control factors have 
been found and analyzed with an ANOVA. In many physical or 
engineering problems two or more variables are related and it is 
of interest to model and explore this relationship. Suppose that 
there is a single dependent variable or response Y that depends 
on k independent or regressor variables, for example x1, x2, ..., 
x3. The relationship between these variables is characterized by 
a mathematical model, called a regression model, that is fit to a 
set of sample data. In some instances, the experimenter knows 
the exact form of the true functional relationship between Y and 
Ks ND, ee es Me 

However, in most cases, the true functional relationship is 
unknown and the experimenter can choose an appropriate func- 
tion in order to approximate it. Low-order polynomial models 
are widely used as approximating functions. Thus, a regression 
model can be written in the following form: 


Y = Bot Bixi + 2x2 t+... + PkXk te (1) 


is called a multiple linear regression model with k regressor vari- 
ables. The parameters £; (j=0, 1, . . ., k) are called the regression 
coefficients. The parameter £; represents the expected change in 
response Y per unit change in x; when all the remaining indepen- 
dent variables x; (ij) are held constant. Thus, the regression 
coefficient of a regressor variable is the sensitivity parameter 
of the dependent variable to changes in the regressor variable 
(independent variable). 

The method of least squares is typically used to estimate the 
regression coefficients in a multiple linear regression model. If 
an experimental session is designed by applying the method- 
ology of the Design of Experiments then the investigated 
dependent variables can be described by regression models. 
Different forms of regression models can be obtained from an 
experimental session according to the chosen design of exper- 
iments, for example first (Eq. (1)) or second-order regression 
models: 


k 


k 
Y = Bo+ > Bixi+ >) Y Bijxixy +S Bjt te (2) 
j=l 


i<j j=! 


Because of the reduced experimental domain of this analysis 
we will deal with results obtained by a first-order design. 
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The data have been first analyzed through the Yate’s method 
(analysis of the factor significance on the dependent vari- 
able); after, first and second-order regression models have been 
obtained. Particular attention has been addressed in the descrip- 
tion of the operation in terms of sensitivity maps. The analysis 
is described in the following paragraphs. 


3. Results 


3.1. Sensitivity of the voltage distribution to fuel 
consumption and air stoichs 


In this paragraph the results of the sensitivity analysis of local 
voltages to fuel consumption and air stoichs are shown. As said 
in the last paragraph, the coefficients of the regression models 
represent the sensitivity coefficients of the analyzed dependent 
variable to the factor. In this case the coefficients are expressed 
in V/%FC and V/%stoichs, being the FC and Agi, factors the 
global values (inputs of the operator). 

In order to better understand the discussion about sensitivity 
to fuel utilization tests some results are presented about the volt- 
age distribution in the stack and the polarization behavior [18]. 
In Fig. 4 (non-disclosure agreement with Siemens) the voltage 
distribution in the stack sectors is shown for the following oper- 
ating condition ([= 435 A, setpoint temperature Toen = 967 °C, 
FC = 84.25%). In Fig. 5 a short polarization (from /=435 to 
475 A) is shown. 

The average cell voltage is around 666+ 14mV, and 
the estimated average area specific resistance is around 
1.04 + 0.08 2 cm?, comparable with other works [16,19]. The 
performance homogeneity (obtained by the estimated standard 
deviation) of the single cell voltage is then 2%, which is a very 
good index of the operation; during different experimental ses- 
sions, values of performance homogeneity of the single cell 
voltage have been found up to 5%. It is also evident that the 
edges’ sectors show a lower terminal voltage. 

In Fig. 6 the sensitivity of the voltage of a sector (group of 
16 cells in electrical series) to fuel consumption is shown. The 
average sector sensitivity is —44 mVsector/%FC. The increase 
of the fuel consumption leads to a significant decrease of volt- 


Voltage distribution in the main stack sectors 
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Fig. 4. Voltage distribution in the main stack’s sectors. 
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Fig. 5. Polarization behavior of two stack’s sectors. 


age in all the stack’s sectors, mainly because of the reduction of 
the average Nernst potential and the increase of diffusion over- 
voltages. The shape of the sensitivity voltage distribution in the 
bundle rows is parabolic, with a minimum in the central zone of 
the generator: the sectors placed at the edges have very high sen- 
sitivity to FC, in particular the sectors placed at the generator’ 
corners of the ejectors’ side. 

In Fig. 7 the sensitivity of the sector voltage to the air 
stoichs is shown. The average sector sensitivity is around 
—20mVsector/0.1 stoichs. The increase of the air stoichiom- 
etry, Agir leads to a significant decrease of voltage in all the 
sectors (except for the ones placed at the generator’s edges: but 
we have evaluated that the effect is statistically not significant). 
Note that the central zone of the generator shows a maximum 
of sensitivity of voltage to air stoichiometry: it is interesting to 
notice the parabolic shape of the sensitivity distribution along 
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Fig. 6. Sensitivity of the sectors’ voltage to fuel consumption. 
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Voltage sensitivity to air stoichiometry 
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Fig. 7. Sensitivity of the sectors’ voltage to air stoichs. 


the bundle rows and the complementary behavior between the 
air stoichs and the FC sensitivity. 

As the SOFC CHP- 100 is an air-cooled system, the decrease 
of voltage, due to an increase of air stoichiometry, is mainly due 
to the variation of the equilibrium temperature: an increase of 
Aair Causes a reduction of the equilibrium temperature with con- 


Voltage sensitivity of the cells placed 
at the power lead side 
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sequent increase of overvoltages (mainly ohmic contribution). 
This decrease of the temperature is almost uniform in all the bun- 
dle rows of the generator (the discussion is addressed in the next 
paragraph), then also at the edges. Thus, it would be expected 
that also the edge sectors will present the same behavior, that 
is a reduction of voltage with increasing of Air. The negligible 
effect of this operation on these sectors is explained with the 
hypothesis that it is the local fuel consumption the dominant 
factor of the operation in the edge sectors. 

The last consideration concerns the linear dependence 
between voltage and temperature expressed by the obtained 
regression models. Through the statistical analysis it was evalu- 
ated that the first-order regression had a good fit (reduce lack of 
fit value). But from a physical point of view, the linear behav- 
ior is not expected and the performance is expected to be an 
exponential dependence of operating temperature. In fact, the 
solid oxide fuel cell behavior can be estimated by the Arrhenius 
expression as this is the behavior of ionic and electronic conduc- 
tivities of the cell layers [20]. Therefore, the linear dependence 
is explained with the small variation of the experimental domain 
for temperature (953°C at air = 4.6 and 944 °C at Agi, = 4.8). 

Looking at Fig. 8, it is interesting to notice the distribution 
along the bundle rows of the voltage sensitivity of the single 
cell to both fuel consumption and air stoichiometry, outlined for 
cells placed at the ejectors and power leads’ sides. 

The estimation of the single cell sensitivity to fuel consump- 
tion is around —2.8 mVcell/%FC, comparable with other works 


Voltage sensitivity of the cells placed 


at the ejector side 
T T T T 


T T T T T T 


1 
=à 
a 


FC Sensitivity mVcell%FC 
à Sensitivity mVcell/0.1 stoich 
DN i 
a N 


-3 


FC 


stoichs 


Fig. 8. Voltage sensitivity to the factors of the single cells along the bundle rows. 
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Fig. 9. Temperature distribution along the bundle rows. 


[13]; the homogeneity of the FC sensitivity among the various 
cells is around 10%. The single cell sensitivity to air stoichiom- 
etry has been estimated in around — 1.3 mVcell/%0.1 stoichs. 
Also from this figure it is clear that the higher is the FC sensitiv- 
ity the lower is the sensitivity to the air stoichiometry. Moreover, 
the distribution of FC sensitivities is used in the following of the 
paper in order to estimate the local fuel utilization of the cell 
sectors. 


3.2. Sensitivity of the temperature distribution and profiles 
to fuel consumption and air stoichs 


In this paragraph the results of the sensitivity analysis of 
temperature distribution to fuel consumption and air stoichs are 
shown. The coefficients of the regression models, that represent 
the sensitivity coefficients of the analyzed dependent variable to 
the factor, are expressed in °C/%FC and °C/%stoichs. 

In Fig. 9 (non-disclosure agreement) the distribution of 
temperatures in the stack is shown, taking into account two mea- 
surement planes: the bottom plane (corresponding to the closed 
end of tubular cell and to the fuel inlet) and the middle plane. 
With respect to the middle plane temperature, the highest tem- 
peratures are measured in the central zone and at the border 
of the stack with a cooling effect at the fuel ejectors’ side; the 
temperature gradient is more than 100°C [21]. 

In order to better understand the discussion about the sensi- 
tivity analysis, some results about the axial shape of temperature 
along the tubular cell are also presented. In Fig. 10 (non- 
disclosure agreement), the typical axial temperature profile 
is shown (at Tajp=612°C and Agir=4.6, ie =200mAcm~?, 


FC = 84.25%). The average top plane temperature is around 
860 °C, the average middle plane temperature is around 930°C 
and the bottom plane temperature is around 840°C. The low 
temperature at the closed end of the cell (corresponding of 
the fuel inlet condition) is explained by the residual reforming 
endothermic reaction which takes place at the cell anode surface 
[22]. 

In Fig. 11 the sensitivity of the local temperature to air sto- 
ichiometry is shown for the two investigated measurements 
planes (bottom and middle plane). The regulation of the air stoi- 
chiometry affects the thermal equilibrium of cells; an increase of 
Aair leads to a reduction of temperature and the effect is stronger 
at the middle plane (in the figure the results are presented in 
absolute value). Moreover, the central zone of the stack is greatly 
affected by this regulation, due to the design of the air distribu- 
tion plenum. Nevertheless, we observe that the sensitivity of the 
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Fig. 10. Axial temperature profile of a tubular SOFC. 
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Fig. 11. Sensitivity of local temperature to air stoichiometry. 


local temperature is almost uniform along the bundle rows: this 
point will be discussed later on. 

In Fig. 12 the sensitivity of the local temperature to fuel 
consumption is shown. The effect is negative (in figure the 
absolute values are reported): it is interesting to notice how 
the middle plane temperature is slightly affected by the oper- 
ation, while the bottom plane temperatures seem to be affected 
by variations on FC values. The test has been performed at 
fixed setpoint temperature. This means that the increase of FC 
would lead to a stack heating (due to irreversibility), but due 
to the imposed setpoint temperature, an higher air mass flow is 
sent to the stack with the final effect that middle plane tem- 
peratures do not change (in the middle plane is located the 
thermocouple which measures the setpoint temperature) while 
a decrease of temperature is measured in the bottom zone of 
cells. 


4. Deduction of the local fuel utilization distribution 


Experiments about sensitivity to fuel utilization form an 
important feature of testing solid oxide fuel cell systems. In this 
paragraph we describe the sensitivity to fuel utilization tests per- 
formed on the SOFC CHP-100 stack and its relation to the local 
fuel utilization of the cell sectors. In fact, the local fuel utiliza- 
tion of cell sectors can be estimated coupling experimentally 
performed sensitivity to fuel utilization tests with a proposed 
model of cell voltage sensitivity to fuel utilization. 

The following results were obtained from the experimen- 
tal analysis: stack terminal voltage sensitivity — 1.06 V/%FC; 


cell sectors sensitivity —44 mVsector/%FC; single cell sensi- 
tivity —2.8 mVcell/%FC. It has been shown in Fig. 8 that the 
homogeneity of the FC sensitivity among the various cells is 
around 10%. The different behavior of the sensitivity is mainly 
addressed to the local fuel utilization operation, so to the design 
of the fuel distribution system. In order to deduce the local fuel 
utilization a model of cell voltage sensitivity to FU has been 
developed. 

The sensitivity of cell voltage to fuel utilization depends on 
several contributions which concern the Nernstian term, the con- 
tribution of diffusion term and eventually the effect of leakages 
of air at the anode side. It is neglected the ohmic contribution, 
the activation contribution is also neglected because of its small 
effect of the voltage drop. If the performance of a cell is limited 
by the ohmic contribution and leakages can be neglected then 
the sensitivity of cell terminal voltage to FU is well described 
by the variation of Nernstian term with fuel utilization. In the 
case of a tubular cathode supported cell the performance is also 
limited by cathodic diffusion and leakages play an important 
rule [13]. In the following some considerations are made about 
all these terms. 

The terminal voltage of a tubular SOFC is expressed by the 

equation: 
Ve = Vnerst — Nact,a/e — Nohm — "Nconc,a/c (3) 
where Ve is the terminal voltage, VNernst the Nernst voltage, 
Nact,a/c the activation overpotential, ohm the ohmic overpoten- 
tial, Nconc,a/c iS the concentration overpotential. 
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Fig. 12. Sensitivity of local temperature to fuel consumption. 


The Nernst voltage is expressed, with the hypothesis of ideal 
gas, by the following equation: 


RT Omo) 
VNernst = Eo + In ( 2/0 (4) 


2F n0) 


Under the assumption that the oxygen utilization is negli- 
gible, that is the oxidant is fed to the cell far in excess of the 
stoichiometrically required amount, the Nernst voltage is a func- 
tion of the operating temperature and of the axial location along 
the cell length, at fixed pressure: 


RT ( Omo) 
Vi ernst 4, T) = Eo(T li 
Nernst(x, T) o(T) + oF o( aT ES 


(5) 


The average Nernst potential is the position averaged Nernst 
potential, that is: 


- 1 ft 
VNernst = = | Vnernst(x) dx (6) 
L Jo 


Basing on the papers [13,22], it is shown that the sensitivity 
of voltage to FU is depending on the average current density due 
to the variation of the consumption of reactants along the cell 
axis with current: the sensitivity decreases with the increase of 
current. In order to model this behavior some assumptions are 
now made. At high current densities it can be assumed that the 
Nernst voltage is linear with the cell axial position: 


(yu, Do)! 
(yH,0(x)) 


RT 
Vnernst(x, T) = Eo(T) + ZF In ( =Ax+B 


(7) 


This refers to an operation with exponential variation of the 
reactants’ partial pressure along the cell length according to the 
equation: 
yi(x) = y;(0) e((1/L) log(yi(L)/yi(O)))x (8) 

The boundary conditions are assumed by knowing the fuel 
composition at the cell inlet, in terms of hydrogen and water 
vapour molar fractions, and the fuel utilization, this allows one 
to find the A and B coefficients. 

At x=0 (closed end), that is at the cell inlet, it can be written: 

yu, (0) = yh, 
yH,0(0) = yho 


(9) 


in the fuel mixture (anode inlet composition). 
At x=L (open end), that is at the cell outlet it can be written: 


yu, (L) = yy, (1 — FU) ae 
yH,o(L) = Yoyo + yh, FU 


The position averaged Nernst voltage at high current densities 
is written in the form: 


(yu, (1 — FU)) vol) 


7 RT 
Wnemst,c(T) = Eo(T) + — In 
Nernst, HC(T) o(T) 4F ( (ymo F Yup FU) 


(1) 


At low current densities, a significant length of the cell would 
be operating close to the exit Nernst potential, that is most of 
the fuel will be consumed very close to the fuel inlet, which is 
the closed end of the cell. The Nernst voltage would decrease 
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sharply from the inlet value to the exit value to the exit value 
within a short length from the inlet. At low current densities, the 
average Nernst potential can be approximated to that of the exit 
Nernst potential, then constant along the cell tube: 


RT (we = aes) 


Nernst,LC(7) o(T) a (mo + ym FU) 


(12) 


Neglecting the diffusion term and the leakages the sensitiv- 
ity of the cell voltage to the fuel utilization can be assumed to 
be equal to the sensitivity of the Nernst voltage to the same 
parameter. With this assumption it is possible to write: 


dV. = dVNernst 
dFU  dFU 


With this assumption the derivative of the terminal voltage 
with respect to the fuel utilization, equal to the sensitivity to the 
fuel utilization, at low current densities is twice the one at high 
current densities, in fact at high current density the sensitivity is 
given by the Eq. (14), while at low current density it is given by 
the Eq. (15). 


(13) 


OV Nernst RT 1 
= J a (14) 
OFU luc 4F |A — FU)(Yp,o + yy, FU) 
3 VNernst RT 1 
= 5 7 (15) 
OFU lic 2F | (1 — FU)OR,0 + y4, FU) 


If the oxygen utilization would be taken into account then 
the Eqs. (14) and (15) will be modified according to (i.e. on Eq. 
(15)): 


IVN RT | 1 


FU  2F | (1 — FUO% o + YA, FU) 


1 (A/y8,)— å 
2 (a — yf, FUXA/y8,) — Ip FU) 


(16) 


Because of the partial pressure of the oxygen at the cathode 
exit can be written as [16]: 


_ À — FUyĝ, 
(A/y0,) — FUyh, 


In the case of the tested air stoichiometry values, à 4.6/4.8, 
this effect can be neglected on the sensitivity curve. This rela- 
tion involves a definition of air excess different than the one 
used in the paper (such values are higher than ones tested exper- 
imentally, see ref. [16]), thus all considerations discussed in the 
paper represent a conservative situation and the discussion is 
only deepen for a comprehensive analysis of all polarization 
terms. 

To refine the model some considerations about the effect 
of fuel utilization on the diffusion losses are now discussed. 
The analyzed cell is a cathode-supported cell then the main 
effect of diffusion is mainly addressed to the cathodic layer. 
The model equations and hypothesis are assumed by literature 


YO» (17) 


[22,23]. Transport of gaseous species usually occurs by binary 
diffusion, where the effective binary diffusivity is a function of 
the fundamental binary diffusivity Dy,—Dy,0 and microstruc- 
tural parameters of the anode. In electrode microstructures with 
very small pore size, the possible effects of Knudsen diffu- 
sion, adsorption/desorption and surface diffusion may also be 
present. In terms of physical measurable parameters, an analyti- 
cal expression for anodic concentration polarization is proposed. 
This equation is: 


2 Fp, (FU)D 
pom i Ae (18) 
RTta 


It is obtained by posing to zero the following equation: 


RTta 


19 
2FDxet (e) 


Where the pressure of fuel at the bulk depends on fuel uti- 
lization operation. 
The voltage drop due to the anodic diffusion is given by: 


Vit = aig 1 (20) 
diff — oF 8 ins 
In terms of physical measurable parameters, the cathode- 
limiting current density can be evaluated in the form: 
— 4F pe, Devett) 
cs — 
(Cp — pQ,)/P)RTte 


This equation is an approximation for anode supported cells 
with very small cathode thickness. In general the limiting current 
at the cathode is given by posing to zero the equation: 


(21) 


Poy = P — (P = PO, (A, FU)) RMP Petey Ne (22) 


Where the partial pressure of the oxidant at the bulk depends 
on fuel utilization and air stoichiometry (Eq. (17)). 
The voltage drop due to the cathodic diffusion is given by: 


> RT i 


In this way the drop of voltage due to a variation of fuel utiliza- 
tion is modeled also considering the effect on the cathodic diffu- 
sion. The effect of diffusion on the terminal voltage sensitivity 
of fuel utilization has significant effect when the limiting current 
density is approached or for low value of the air stoichiometry. 

It is possible to write with respect to Eqs. (20) and (23): 


Vie _ RT 1 (nN 


əFU  2F1 — (ie/ias) A 
2F Cias = ic)las dFU 
OVS RT 1 ic(Oics/JFU) 
ðU 4F1— (iofice) i2 
_ RT ic s (25) 


AF (ics — ig)ics OFU 
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Table 1 
Parameters used in the modeling of the diffusion overvoltages 


Parameter 


Anode 


Daurceft), effective binary diffusivity 


D; — Dj, fundamental binary diffusivity (cm? s7!) 


T 3 


293 K = 
DK o S051 


E 0.3 


Cathode 
Daet = (€/t) Dip, $p o(T/293)!5 (e/D DRX (1/293)! 
293K 
DO, Ny = 0.22 
0.335 
5 


It should be noted that: 


=a/C 
avait l 


—— (26) 
dFU Da/c(etf) 
ayas 
r T X fajo (27) 


The effect of the diffusion on the sensitivity of voltage to the 
fuel utilization is inversely proportional to the effective binary 
diffusivity on the electrodes’ layers and directly proportional to 
the electrodes’ thickness. 

The parameters used for solving the equations have been 
listed in Table 1. 

An important role is played by eventual leakages of air at 
the anode side; in fact some fuel can burn with the oxidant with 
the consequence that it does not react at the anode surface to 
generate current. If it would occur then the effective operational 
fuel utilization will be higher than the expected and sensitivity 
to change in FU will be very high approaching high values of 
fuel utilization. 

The results of modeling are shown in Fig. 13. The simplified 
model refers to the sensitivity model of Nernst voltage to fuel 
utilization whereas the complete model accounts also for the 
effect of diffusion and leakages (3% of leakages). 

The effect of the operating current on the fuel utilization sen- 
sitivity is more effective when working at low fuel utilization. At 
low current the fuel utilization sensitivity is governed essentially 
by the variation of the average Nernst voltage with FU whereas 
at high current the effect of diffusion and leakages becomes 
more effective. The curves are drawn for an operating condi- 
tion of 1000°C and for a fuel with 89% Hz and 11% H20 of 
composition. In Fig. 14 some experimental results of single cell 
sensitivity are reported according to [13]. 


Model of voltage sensitivity to FU 
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Fig. 13. Model of fuel utilization sensitivity for a tubular SOFC. 
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Fig. 14. Experimental data of single cell sensitivity. 


Once the model is developed, it can be used in order to explain 
the different values of local sensitivity to fuel utilization and 
to analyse the shape of its distribution along the stack sectors 
(estimated experimentally through the statistical analysis). The 
distribution of fuel utilization sensitivity (according to the pro- 
posed model) depends on the local temperature in the bundle 
rows, on the fuel inlet condition at the anode in terms of hydrogen 
and water vapour partial pressures, and on the local fuel utiliza- 
tion. The effect of the local temperature is taken into account in 
the analysis, the effect of the fuel inlet conditions is taken into 
account but the strong hypothesis is that the fuel composition is 
the same for all the bundle rows, that is there is the same ref- 
ormation degree in all the in-stack standard reformation boards. 
Once the temperature and the fuel inlet condition are fixed for all 
the cell sectors, then the local fuel utilization can be estimated by 
solving the voltage sensitivity model in terms of fuel utilization. 
The condition of the inlet fuel has been posed to yp, = 0.4 and 
ymo = 0.2. The results of the estimation procedure are shown 
in Fig. 15; the higher is the local FU sensitivity the higher is 


Estimation of the local fuel utilization 
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Fig. 15. Estimation of local fuel utilization. 
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g. 16. Model of FU voltage sensitivity vs. estimated local fuel utilization. 


the local fuel utilization, the boundaries of the generator work at 
high fuel utilization, and this can also explain the lower voltages 
observed at the edge’s sectors. 

In Fig. 16 the deduced local fuel utilization with the FU 
sensitivity model used for the estimation are compared. 

The model expresses the FU voltage sensitivity at the nominal 
current of ie =200 mA cm~?. It is shown that the difference of 
local fuel utilization of cell sectors can differ also for 10% of 
values. 


5. Air stoichs regulation to avoid local overheating: 
some problems 


The model discussed in Section 4 has been used to under- 
stand why there is a distribution of voltage sensitivities to FC 
in the stack as shown in Fig. 6, and to deduce the local FU 
distribution in the various sectors of the stack: the higher is the 
sensitivity to fuel utilization the higher is the estimated local fuel 
utilization. The distribution of the local FU value in the 12 rows 
of the generator is therefore shown in Fig. 17 (non-disclosure 
agreement). 

The distribution shows higher values of local FU in the bound- 
ary rows of the generator, and lowest (better) values in the center 
rows. This is probably due to the fuel gas distribution in the gen- 


Relation between local temperature and fuel utilization 
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Fig. 17. Distribution of local FU and temperature in the rows of the stack. 


erator rows and to the operation of the boundary reformation 
boards: it seems that in the boundary rows a lower amount of 
fuel gas arrives and the operation is characterized by a local FU 
factor higher than the overall one. For this reason lower voltages 
are also observed in these regions (see Fig. 4). 

In the same Fig. 17 the local temperature distribution is 
reported. We observe that it seems to exist a relation between the 
local FU and the local temperature of the rows: the higher is the 
local FU the higher is the local temperature. Moreover, it could 
be interesting to consider the fact that the boundary rows (1 and 
12) have the reformer unit on the external side characterized by 
an operation with a reduced amount of fuel gas, compared to the 
reformation boards placed internally between the inner bundle 
rows. The amount of the reforming reaction is therefore reduced 
in these boundary reformers, and therefore its cooling effect is 
reduced. This could cause the higher values of both local FU 
and local temperature. 

Thus, it is outlined how high local fuel utilization values can 
lead to local overheating of the generator. Of course, the local 
overheating causes some problems. First, it limits the opera- 
tion of the generator, especially the value of current density that 
the generator can reach, because an increase of current density 
determines an increase of the local temperature with too high 
values at the boundary rows. The limitation in the current den- 
sity means a limitation in the power density of the generator, 
which is one of the most important parameters of the genera- 
tor performance. Second, the local overheating can accelerate 
the degradation of the materials of the tubular cells, causing the 
degradation of the generator performance and the necessity of a 
increased maintenance. 

To overcome the problem of the local overheating, the first 
attempt is to analyse the behavior of the system devoted to the 
thermal regulation of the generator. The thermal regulation is 
assured by the cathode air flow, which is distributed inside the 
generator rows. The thermal regulation is made by modifying 
two factors: the air stoichs and the air pre-heating temperature, 
but the most important factor is the air stoichs. The air distri- 
bution is made with a unique air plenum, which means that the 
modification of the air stoichs is not differentiate for the differ- 
ent segments of bundle rows, but it is common for the whole 
generator and the various bundle rows. Therefore, through the 
use of the sensitivity analysis of the effect of the air stoichs on 
the local temperature (described in Section 3.2 and shown in 
Fig. 11), we have determined the sensitivity of the local temper- 
ature to a modification of the overall air stoichs factor: the effect 
is shown in Fig. 18. 

The results are shown for the ejectors side, the central zone 
and the power lead side of the generator. The sensitivity is higher 
in the central zone (as reported in Fig. 11), but the most important 
notation is that the sensitivity is almost constant along the 12 
bundle rows of the generator, that is, it does not depends on the 
position of the bundle row inside the generator. Therefore, we 
can expect that a modification of the air stoichs determines a 
similar effect in all the bundle rows, that is, the cooling effect 
is not differentiate in the different zones of the generator. In 
fact, a modification of the temperature distribution consequent 
to a regulation of the air stoichs (increase of the air stoichs from 
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Fig. 18. Sensitivity of the local temperature to a air stoichs regulation. 


Aair = 4.6 to Agir = 4.8, that is, a cooling of the generator) is shown 
in Fig. 19 (non-disclosure agreement). 

As we see, the thermal regulation through the increase of the 
air stoichs determines a reduction of local temperature almost 
uniform along the bundle rows. The regulation does not solve the 
problem of the local overheating at the boundary rows, and has 
the consequence of a uniform reduction of temperature, with a 
decrease of the mean operation temperature of the generator, and 
a consequent increase of the overvoltages. Therefore, the present 
air cooling system is able to change the thermal level of the stack, 
but in a uniform way and cannot act at the local level; therefore it 
cannot control a local overheating. This limits the possibility of 
increasing the setpoint temperature, to reduce the overvoltages. 
But the main drawback is the difficulty of increasing the current 
density value, with a consequent limitation of the power density 
of the generator. 

An expected thermal regulation of the generator is shown 
in the figure: a regulation able to obtain an expected uniform 
value of temperature equal to the setpoint value imposed by the 
operator. This regulation could be able to overcome the local 
overheating, while assuring that the mean temperature main- 
tains its imposed setpoint value. This regulation can be obtained 
through a diversification of the air stoichs imposed in the differ- 
ent zones (or in the different bundle rows) of the generator. The 
modification has to be made at the level of the air distribution 
systems, with different air plenums which could differentiate the 
air stoichs regulation for the different bundle rows. 


Effect of the temperature regulation by the present air cooling system 
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Fig. 19. Effect of the temperature regulation by the present air cooling system. 


6. Conclusions 


In the paper, an analysis of a SOFC system is presented in 
terms of sensitivity maps of local voltages and temperatures with 
respect to the overall fuel consumption and the air stoichiometry 
control factors: 


e the local voltage sensitivity tests to fuel consumption are used 
as a useful procedure to estimate the local fuel utilization of 
the main stack’s cell sectors during the operation: the local 
fuel utilization distribution has been found; 

e high local fuel utilization have been detected: they can lead 
to local overheating which can limit the operation (power 
density, efficiency and degradation) of the generator; 

e with the present air cooling system design it is possible to vary 
uniformly the operating temperature of the generator, but it is 
not possible to control the local temperature; 

e the present air cooling system is able to change the thermal 
level of the stack, but in a uniform way and cannot act at 
the local level; therefore it cannot control local heating: this 
limits the possibility of increasing the current density with a 
consequent limitation of the power density of the generator; 

e amore effective regulation should be able to obtain a uniform 
value of temperature equal to the setpoint value imposed by 
the operator; this regulation can be obtained through a diver- 
sification of the air stoichs imposed in the different bundle 
rows of the generator, through the adoption of different air 
plenums which could differentiate the air stoichs regulation 
for the different bundle rows. 


Acknowledgement 


This work was financially supported by Gas Turbine Tech- 
nologies (Contratto Attuativo 2006), in the frame of the EOS 
Project. 


References 


[1] S. Saroglia, P. Leone, M. Santarelli, Eco-efficiency 2005-HYSYDays Ist 
World Congress of Young Scientists on Hydrogen, Torino, Italy, May 
18-20, 2005. 

[2] P. Leone, M. Santarelli, M. Cali, J. Power Sources 156 (2006) 400-413. 

[3] M. Cali, M.G. Santarelli, P. Leone. Int. J. Hydrogen Energy, in press. 

[4] P. Leone, M. Santarelli, Eco-efficiency 2005 — HYSYDays Ist World 
Congress of Young Scientists on Hydrogen, Torino, Italy, May 18-20, 2005. 

[5] M. Cali, M. Santarelli, P. Leone, Comparison of the Behavior of the CHP- 
100 SOFC Field Unit Fed by Natural gas or Hydrogen through a Computer 
Experimental Analysis, World Hydrogen Technology Convention, Singa- 
pore, 2005. 

[6] M. Cali, E. Fontana, V. Giaretto, G. Orsello, M. Santarelli, HySafe— 
International Conference on Hydrogen Safety, Pisa, Italy, 2005. 

[7] D.C. Montgomery, Design and Analysis of Experiments, John Wiley & 
Sons, Inc., 2005. 

[8] R.E. Walpole, R.H. Myers, Probability and Statistics for Engineers and 
Scientists, Prentice Hall International Inc., New Jersey, 1993. 

[9] M. Cali, F. Debenedictis, P. Leone, M. Santarelli, V. Verda, ASME Ist 
European Fuel Cell Technology and Applications Conference, Rome, Italy, 
December, 2005. 

[10] M. Cali, P. Leone, M. Santarelli, G. Orsello, G. Disegna, General 
Topics and Operation Description by Means of Regression Models, 


168 


11 


M. Santarelli et al. / Journal of Power Sources 171 (2007) 155-168 


World Hydrogen Energy Conference WHEC 2006, Lyon, France, June, 
2006. 

M. Cali, G. Orsello, M. Santarelli, P. Leone, Proceedings of ESDA2006, 8th 
Biennial ASME Conference on Engineering Systems Design and Analysis, 
Torino, Italy, July, 2006. 

M. Santarelli, P. Leone, M. Cali, F. De Benedictis, Lucerne Fuel Cell Forum 
2006, Lucerne, Switzerland, July, 2006. 

S. Gopalan, G. DiGiuseppe, J. Power Sources 125 (2004) 183-188. 

M. Molinelli, D. Larrain, N. Autissier, R. Ihringer, J. Sfeir, N. Badel, O. 
Bucheli, J. Van Herle, J. Power Sources 154 (2006) 394-403. 

R.A. George, J. Power Sources 73 (1998) 251-256. 

S.C. Singhal, K. Kendall, High Temperature Solid Oxide Fuel Cells: Fun- 
damentals, Design and Applications, Elsevier, 2004. 

H. Kabs, Operational Experience with Siemens-Westinghouse SOFC 
Cogeneration Systems, Lucerne Fuel Cell Forum, 2001. 


[18] M. Santarelli, P. Leone, M. Gariglio, M. Cali, P. Spinelli, G. Orsello, F. De 
Benedictis, Fuel Cell Seminar 2006, Honolulu, Hawaii, USA, 2006. 

[19] S.C. Singhal, in: W. Poulsen, N. Bonanos, S. Linderoth, M. Mogensen, et al. 
(Eds.), Proceedings of the 17th Riso International Symposium on Materials 
Science: High Temperature Electrochemistry: Ceramics and Metals, Riso 
National Laboratory, Roskilde, Denmark, 1996, p. 123. 

[20] K.S Huang, Cell Power Enhancement via Materials Selection, 7th European 
SOFC Forum, Lucerne, Switzerland, July 3-7, 2006. 

[21] S.D. Vora, Advances in Solid Oxide Fuel Cell Technology (SOFC) at 
Siemens Westinghouse, Fuel Cell Seminar, Palm Springs, California, USA, 
November 14-18, 2005. 

[22] R. Draper, G. DiGiuseppe, Proceedings of ASME European Fuel Cell 
Conference, Rome, Italy, December 14-16, 2005. 

[23] J. Kim, A.V. Virkar, K.Z. Fung, K. Metha, S.C. Singhal, J. Electrochem. 
Soc. 146 (1) (1999) 69-78. 


